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We use a phenomenological Hamiltonian approach to derive a set of coupled mode equations that 
describe light propagation in waveguides that are periodically side-coupled to microcavities. The 
structure exhibits both Bragg gap and (polariton like) resonator gap in the dispersion relation. The 
origin and physical significance of the two types of gaps are discussed. The coupled-mode equations 
derived from the effective field formalism are valid deep within the Bragg gaps and resonator gaps. 
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I. INTRODUCTION 

In the past several years the linear and nonlinear prop- 
erties of side-coupled waveguiding structures have at- 
tracted the attention of many researchers 1 - - — . These 
structures consist of one or more waveguiding elements 
in which forward and backward propagating waves are 
indirectly coupled to each other via one or more mediat- 
ing resonant cavities. Perhaps the most common propos- 
als for realizing these structures involve photonic crys- 
tal (PC) waveguides with defect modes slightly displaced 
from the waveguiding region (Fig. la, left)£i&, or micro- 
ring resonator structures in which two channel waveg- 
uides are side-coupled to micro-ring resonators (Fig. la, 
right) ». In the PC structure the forward and backward 
propagating modes within the waveguide are coupled via 
the defect; for the micro-ring structure, the forward going 
mode in the lower (upper) channel waveguide is coupled, 
via the micro-ring, to the backward going mode in the 
upper (lower) channel. The linear and nonlinear proper- 
ties of both types of structures have been studied 2 ^ - — . 

The electromagnetic properties of these structures can 
be accurately determined in great detail using numer- 
ically intensive methods such as finite-difference time- 
domain (FDTD) simulations 14 . An analysis in terms 
of Wannier functions can substantially reduce compu- 
tation time for the PC structure 1 ^, but the numerical 
problem remains daunting. In particular, full FDTD cal- 
culations of the micro-ring structures have to date been 
confined to two-dimensional analogs of the actual struc- 
tures of interest* 1 ^. Furthermore, direct numerical sim- 
ulation, while valuable for design purposes, offers little 
insight into the physics of the structures. Consequently, 
semi- analytical techniques, such as the scattering-matrix 
approach of S. Fan et alrk and Yong Xu et al&, have been 
proposed. Using these techniques the optical properties 
of side-coupled structures can be understood in terms of 
the interactions between a small number of modes. 

In this paper we concentrate our attention on periodic, 
side coupled structures (Fig. lb). Our primary objective 



is to derive coupled mode equations (CME) that describe 
pulse propagation in such structures. Coupled mode the- 
ory has long been used as an effective design tool for grat- 
ing structures where forward and backward propagating 
waves are directly coupled via an index grating— In di- 
rectly coupled structures, it is well known that a Bragg 
gap opens in the dispersion relation of the structure when 
the phase accumulated in one round trip through a pe- 
riod of the grating is an integer multiple of 27r, so that 
the slight reflections that are incurred due to the grat- 
ing are coherently enhanced. Structures possessing a 
Bragg gap have found a variety of uses, such as dispersion 
compensation^ and wavelength division multiplexing 18 . 
In the side-coupled structure the Bragg feedback mech- 
anism, and hence the Bragg gap, does exist, although it 
is now mediated by the coupling cavity. However, there 
is also a second type of gap: a resonator gap, which is 
associated with the resonance frequencies - and therefore 
the geometry - of the mediating cavity. For the micro- 
ring resonator structure the interpretation of this gap is 
straightforward: when the phase accumulated in a round- 
trip through the micro-ring resonator is an integer mul- 
tiple of 2n, then the coupling between the forward and 
backward going waves is resonantly enhanced. Of these 
two gaps, the resonator gap is perhaps the more impor- 
tant, because it exhibits a deep transmission dip seen 
even in a structure with only one unit cell. 

Because side-coupled structures exhibit both Bragg 
and resonator gaps, it is to be expected that a CME de- 
scription of optical pulse propagation will be more com- 
plicated than in Bragg gratings. The CME for Bragg 
gratings involve two fields (forward and backward going) 
interacting via a coupling coefficient. For side-coupled 
structures, the most interesting situation is when a res- 
onator gap lies near one of the Bragg gaps, and we show 
in this paper that the relevant CME then involves three 
fields: a cavity field and forward and backward going 
fields. 

We derive our CME using a phenomenological Hamil- 
tonian approach, which distills the essential physical in- 
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teractions of the structure, and hence provides a simple 
physical picture of optical interactions. We build the 
fields in our CME as Fourier superpositions of the modes 
in the Hamiltonian. Hence, our CME are derived for in- 
finite, periodic structures in which the coupling to each 
cavity is the same. Nevertheless, we show that our CME 
can be generalized to describe finite, apodized structures, 
in which the coupling (but not the period) varies from 
cavity to cavity. Therefore, the CME can be used to 
describe finite structures with only a small number of 
cavities. Indeed, the general Hamiltonian approach we 
advocate can be applied even to structures with only one 
or two cavities, if the formalism we introduce in Sec. II 
is extended to a discrete number of (not necessarily iden- 
tical) cavities. In both discrete and periodic scenarios, 
the Hamiltonian approach exhibits the similarities of the 
optical dynamics of these artificially structured materials 
to more traditional problems in solid state physics. As 
well, it allows for an easy quantization of the description 
to address the quantum optics of these structures. We 
plan to turn to this, as well as the direct derivation of 
our phenomenological Hamiltonian from the underlying 
electrodynamics, in future publications. 

The present paper is organized as follows. In Sec. II we 
describe the Hamiltonian model for a system with a single 
microresonator, investigate the transmission/reflection 
spectrum of the structure, and indicate how the param- 
eters in our phenomenological Hamiltonian can be set 
from more common models of cavity resonators. In Sec. 
Ill we discuss how the Hamiltonian can be used to model 
a periodic waveguide-resonator structure. We then dis- 
cuss methods of reducing the number of fields and in- 
teractions in our Hamiltonian while retaining the basic 
physics. In Sec. IV we derive the coupled mode equations 
in terms of effective fields built as Fourier superpositions 
of the modes in the Hamiltonian of Sec. Ill, and we show 
how to modify these CME to describe finite, apodized 
structures. In Sec. V we conclude. 



II. HAMILTONIAN MODEL AND 
TRANSMISSION FOR A SINGLE CAVITY 
STRUCTURE 

In this section we construct a Hamiltonian model for 
a structure in which forward and backward propagating 
waves are indirectly coupled to each other via a cavity 
centred at z = zq. We will focus on classical optics here, 
but because its easy generalization to quantum optics is 
one of the strengths of this approach, we adopt a quan- 
tum notation and, for the classical Poisson bracket {.., ..}, 
we write (ift) -1 [.., ..]; we also use t to indicate complex 
conjugation. We will also often speak of operators rather 
than variables, especially when it makes the physics more 
clear. For example, we introduce a\ and c\ as creation 
operators for photons propagating with wavenumber k 
in the forward and backward direction respectively. Be- 
cause k > (k < 0) indicates that the photons are prop- 
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FIG. 1: (la) Waveguide-resonator structure containing (left) 
photonic crystal microcavity and (right) micro-ring resonator. 
On the right the micro-ring resonator is coupled to two waveg- 
uides, with the forward (backward) propagating light in the 
lower (upper) waveguide; On the left the structure contains 
dielectric rods is embedded in air, the singly degenerate mi- 
crocavity is coupled to the photonic crystal waveguide, formed 
by removing a row of rods in the photonic crystal, (lb) Peri- 
odic waveguide resonator structure containing (top) microring 
resonator and (bottom) photonic crystal microcavity. 



agating in the forward (backward) direction, at exists for 
k > and c\ for k < 0. For a given k, the energy in these 

fields is hoJku k ak and fru>kC k Ck, with u>k — c \ k\ /n, where 
c is the speed of light in a vacuum, and n is a constant 
effective index, equal for the forward and backward prop- 
agating waves. By ignoring the frequency dependence of 
n we are neglecting the underlying material dispersion 
within the waveguides; we discuss the validity of this ap- 
proximation after equation @ below. To describe light 
in the cavity, we define a creation operator , and iden- 
tify the energy in the field as hujQb^b, where luq is the 
resonant frequency of the cavity. For the micro-ring res- 
onator structure of Fig. la (right), the a\ and c k could 
represent creation operators for light propagating in the 
forward direction in the lower waveguide and the back- 
ward direction in the upper waveguide, while could 
represent the field circulating in the counter-clockwise 
direction in the micro-ring resonator. Our notation im- 
plies that the two waveguides have a common mode in- 
dex n, but this could easily be generalized. For the PC 
structure of Fig. la (left), the a k and c k would represent 
creation operators for light propagating in the forward 
and backward direction in a waveguide mode of the PC 
waveguide, and b^ would represent the creation operator 
for the field inside the single mode defect. Regardless of 
their interpretation, the operators satisfy the commuta- 
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tion relations 



a,k,a k , 
Ck,c[, 



= S(k-k') 



= 6(k-k') 
[6,6+] = 1, 



(1) 



with all other commutation relations vanishing. Assum- 
ing that no light couples directly between the propagat- 
ing modes governed by at and ct, but that light can cou- 
ple from these modes to the cavity, we use the following 
model Hamiltonian for the systern^^: 



H=H Q + H, 



coupling 3 



(2) 



where 



dkhujka k ak 



r° 

/ dkfiujkC k Ck + huiob^b, 
J — 00 



H, 



coupling 



-n I & 
10 



a\be- ikzo + 6+a fc e* zo 



dk 



{-l) q h £_Jc+6e- ifez ° + tfc k e ikz ° 



(3) 

dk. 

(4) 



The quantities and (— l) 9 characterize the strength 
of the coupling between cavity field and waveguide fields, 
propagating in the forward and backward direction; q is 
an integer that depends on the symmetry of the cavity 
mode&. Note that except for the factor (— l) 9 our nota- 
tion implies that the coupling to forward and backward 
propagating waveguide modes is identical. In the micro- 
ring structure, for example, this means that we assume 
equal coupling to the two waveguides; generalization of 
this is straightforward, but for simplicity we will not do 
it here. The time evolution of the operators is given by 
the Heisenberg equations of motion 



(5) 



where O is any operator. 

In writing down J5J, © and Q we have implicitly 
assumed that the cavity supports only one mode, with 
resonant frequency loq, and that the waveguides guide 
light in only a single spatial mode profile. Strictly speak- 
ing, of course, neither of these assumptions is valid. In 
general, cavities support more than one mode, oscillating 
at one or more resonance frequencies, and for sufficiently 
high frequencies a waveguide will support multiple trans- 
verse modes. However, we are primarily interested in the 
physics of these structures for frequencies at or near a 
specific resonant frequency wo- We then assume that 
within this frequency range only one resonance of the 
cavity exists or, alternatively, that only a single mode of 
a multi-mode cavity is excited, and that the waveguides 
of the structure are single mode. Furthermore, we as- 
sume that the underlying material or modal dispersion 



of the structure is negligible within the frequency range 
of interest. For our purposes, the inclusion of material 
dispersion would lead to quantitative, but not qualitative 
changes. 

In Appendix 1 we show that our Hamiltonian formu- 
lation leads to a Lorentzian transmission and reflection 
across the cavity for frequencies in the vicinity of loq: 



t(w) 



-iA 



7 — iA 



7 — iA 



(6) 
(7) 



where 7 = 2nn^ 2 g /c, and £o D is the coupling coefficient 
between the cavity and waveguides evaluated at k = uj = 
noJo/c, and where A = (lu — luq — a (to)) characterizes 
the detuning from the renormalized resonance frequency 
ujQ + a (lu) . An expression for the quantity a (lu) is given 
in Appendix 1. For our structures of interest a(uj) is 
sufficiently small that u) — lu$ — a (cj) ~ lu — loq to a good 
approximation. 

The transmission and reflection coefficients in ©,0 
are of precisely the form that follows from simple transfer 
matrix models of resonant cavities or ring resonators&i. 
In the latter structure, for example, the coupling of the 
cavity to the waveguides is described by self-coupling and 
cross-coupling coefficients a and k respectively, which in 
a simple case (where the coupling is assumed to occur 
at the point of smallest separation) are real and satisfy 
a 2 + k 2 — 1. Comparing the transmission and reflection 
coefficients found there with |JBJ|,(J2J|, we find that they 
become equivalent if we put 



2nnR 



1 



(8) 



where n and R are the effective index and radius of the 
resonator respectively. Thus if a given resonator is pa- 
rameterized by a and well of course by the resonance 
frequency ojq, then relation JSJ allows one to determine 
the effective coupling coefficient ^ and thus set what 
will be, as we will see, the crucial elements in the phe- 
nomenological Hamiltonian The appropriate values 
of (7 and k for a single resonator could be determined by 
experiment, or directly calculated from the underlying 
channel and resonator geometries, as discussed by Waks 
and Vuckovioi^. 

A typical spectrum for a single cavity structure is 
shown in Fig. 2. On resonance, the reflection induced by 
the cavity reaches 100% (albeit only for a single wave- 
length), and remains significant as long as the detuning, 
A, is on the order of 7. The width of the spectrum is 
dictated by 7, and the larger the coupling to the cavity, 
the broader the resonance. In physical terms, this means 
that as the waveguides are brought closer to the cavity 
of Fig. la, the resonance width increases. 
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FIG. 2: Transmission (solid line) and reflection (dotted line) 
spectrum for the one cell structure obtained using equations 
© and iJTJ. The structure can demonstrate 100% reflection 
and 0% transmission when the frequency is matched to the 
resonance frequency of the microresonator. For comparison 
with later plots, the frequency is normalized with a distance 
A, which we use as the distance between resonators when we 
consider a periodic array. 



III. HAMILTONIAN FOR A PERIODIC 
STRUCTURE 



We now generalize the single-cavity Hamiltonian to 
describe a periodic structure, in which the forward and 



backward propagating modes are coupled to an infinite 
series of periodically spaced cavities (Fig. lb). We as- 
sume that the resonators are not directly coupled to each 
other, although of course they do couple indirectly via 
the waveguides. Generalizing the Hamiltonian J5J to in- 
clude the periodic sequence of resonators, we write 



H = 



/ dkhoJk<i k ak + / dkhojkc\.Ck + hu} bjb, 
Jo J-oc l 

r°° r i 
- h J2J dk ^ k [ b i a ke tkzi + a{b ie - lkzi \ 

-(-l^hY, I i-kdk \b\c k e lkz < +clb ie - ikz '~ 



(9) 

where a\ {^Cj are again the creation operators for light 
propagating the forward (backward) direction. The main 
difference between (0 and J5J) is that we have now in- 
cluded a countably infinite number of resonators, each 
with the same resonance frequency, loq, and associated 
with the creation operator bj, where / indexes the res- 
onator. The resonators are evenly spaced at zi = /A, 
which gives a fundamental reciprocal lattice vector Go = 
2ir/A. The Hamiltonian can be re- written as 



J 



H = 



/ dkhuj k+G a\ +G a k+G + £ / dk%u)k-Gcl_ G Ck-G + £ ^A* 6 * 

s~i J B.Z. s~t J B.Z. i 



l G 



dk^k+G 



B.Z. 



l G 



B.Z. 



I 



(10) 



where ^2 G represents the summation over an infinite 
number of positive reciprocal lattice vectors (with G — 
0, Go, 2Go, ...), and where in the integrations we restrict 
the wavenumber k to the first Brillouin zone (— Go/2 < 
k < Go/2); We sum only over the positive reciprocal 
lattice vectors so that a\ +G and c\_ G retain their associ- 
ation with forward and backward propagation modes re- 
spectively. The operators satisfy commutation relations 



a k+G ,a k , +G , 

c k - G ,4'-G> 
bub), 



= 5(k-k')S G>G > 
= S(k- k') d G>G > 
= 5u>, 



(11) 



with all other commutators vanishing; the first two of 
these follow immediately from |QJ . Because the system is 



periodic, we can identify a countably infinite set of Bragg 
frequencies in <|10|l . These are the frequencies 0Jk±G eval- 
uated at k = or Go /2. Hence, since uJk±G = c\k ± G\ /n 
for ring resonator structures, the M th Bragg frequency 
occurs at oj^^ = M (cGq jlri) (with M > an integer) . 
To simplify l|l(J|) . we introduce the collective operator 



-ikzi 



(12) 



where k is now a continuous variable that ranges over 
the first Brillouin zone. In Appendix 2 we introduce this 
operator by first considering only excitations of the res- 
onators periodic over a length L = iVA, and then taking 
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TABLE I: Parameters used for dispersion relation calculation. 



Physical parameters 


cr = 0.98 


A = 32.0^m 


n = 3.0 


2-kR = 26.3^m 


Numerical parameters 


f = 0.0023/im- 1 


i£onA _ 124 _ 156 

7TC 


" b " A - 124.0 

7YC 


BnA = q Q7 

7TC 



A" — > oo. We find in that limit 



2J hcu bjbi 



B.Z. 



dkhu) b k bk, 



for fc and k' in the first Brillouin zone, with all other com- 
mutators vanishing. In terms of this collective operator 
the Hamiltonian (|10|l becomes 



and that 



b k ,bl] =S{k-k') 



H = dkhuj k+G al +G a k+G + V" / dkhuj k - G c\_ G c k - G + / 

r, J B.Z. r, JB.Z. JB.Z. 



dkhuj b\b k 



dkE+ k 



k+G 



B.Z. 



b\a k+G + a\ +G b k - ^5"! / 

J n JB.Z. 



k+G 



b\c k - G + c\_ G b k 



(13) 



where ^± k ± G = y j?±t±G- I n Table I we give typical 
values for parameters characterizing side-coupled struc- 
tures, and we use them in our sample calculations below. 
There and for the rest of this paper we assume that the 
coupling S± k ± G t is approximately constant at wavevec- 
tors corresponding to frequencies within our region of 
interest, and take Z±k±G' ~ 3. This approximation is 
reasonable if the G of interest satisfy G <C A/c, where 
Afc is the range over which the varies significantly. We 
can expect Afc w 27r/(lj«m) for the structures of interest 
(see Appendix 1), and since G is at most a few times 
Go = 27r/A (= 27r/(32/im) from Table I), this inequality 
is indeed satisfied. 

The dispersion relation of the system can be deter- 
mined by traditional transfer matrix methods, using JSJl, 
@ for the transmission and reflection coefficients of a sin- 
gle resonator. However, to see the connection with the 
coupled mode equations we will derive, we consider deter- 
mining the dispersion relation directly from the Hamilto- 
nian (|13|) , by applying the Heisenbcrg equation of motion 
to generate equations for the time derivatives of a k+G , 
Ck-G and bk- Assuming harmonic time dependence e~ tut 



for the operators, we determine an expression for to as 
a complicated function of the countably infinite set of 
w ±|fc|±G> and the discrete value loq. Alternately (and 
equivalently) we can exhibit the Hamiltonian in a matrix 
form (pJJ 

h = h [ dk ft • y k ■ f fc , (w) 

JBZ 

where 

ffc = ( a !+G ' a fc+2Go> C fc-G ' C !-2G ' ^fe) ' (-^) 

and Vfc contains all of the interactions between the ak+G-, 
Ck- G and b k . Then, by diagonalizing the (infinite- 
dimensional) matrix \l k we can in principle determine 
the dispersion relation of the structure. In Fig. 3 we 
consider a typical uncoupled (in the limit where S = 0) 
and coupled dispersion relation for the structure. The 
dotted line shows the uncoupled dispersion relation, and 
the solid line shows the dispersion relation of the coupled 
system, as determined by the transfer matrix approach. 



If one of the Bragg frequencies is close to the resonant 
frequency u>o, then we show below that a truncation of 



the matrix V k to three terms is a good approximation. 
The restricted Hamiltonian that results is 



6 



H 



dk 



B.Z. 



"k+G' ^k-G' 



LOk+G' 3 



dk+G' 
C k -G' 

bk 



(16) 



where G' is the reciprocal lattice vector associated with 
the forward (backward) band that has cok+G' i^k-G') 
closest to ujo. Here we have assumed that the reso- 
nant frequency is very close to a Bragg frequency with 
its associated gap at the Brillouin zone centre, and so 
log' = u-g> = ^6, where tub is the Bragg frequency clos- 
est to the resonance frequency 8 . We refer to eqn. 16 
as the "three mode model." Its validity near a resonance 
frequency for any particular structure can be formally in- 
vestigated by including the omitted terms in a multiple 
scales analysis, or by simply comparing the dispersion 
relation following from eqn. 16 with a full solution of 
the dispersion relation using a transfer matrix approach. 
This is done in Fig. 4, using the parameters in Table 
I as was done in Fig. 3. In Fig. 4 we also plot the 
imaginary part of k within the gaps. Note that the ex- 
act solution and that from the three mode model are in 
good agreement for the frequency range shown in Fig. 4. 
Such agreement fails at other Bragg frequencies that are 
further from the resonant gap, of course, since the three 
mode model (eqn. 16) only contains the physics of the 
Bragg gap closest to luq. It is to frequencies near ujq that 
we henceforth restrict ourselves. 



IV. COUPLED-MODE EQUATIONS IN THE 
THREE-MODE MODEL 

In this section we derive a set of coupled-mode equa- 
tions which describe pulse propagation in the periodic 
structure, based on the three-mode Hamiltonian 116fl . 
We then demonstrate that although these coupled mode 
equations are derived for an infinite periodic system with 
equal coupling at each resonator, they can, with only 
slight modifications, be used to describe finite systems 
with varying coupling at each resonator. We start by 
defining effective fields in terms of the amplitudes ak+c , 
c k - G , and b k : 



g- (z,t) = 
b(z,t) = 



dk 

B.Z. \/27r 

dk 



dk+G'C , 



b.z. v2tt 
dk 



Akz 

Ck-G'e , 



B.Z. 



<2tt 



--b k e l 



(17) 



where G indexes the reciprocal lattice vector that is 
retained within the three mode approximation. These 
fields can be interpreted as a forward propagating field, 
a backward propagating field, and the field distribution 
in the resonators respectively. Using the definitions in 
p7[l. the effective fields satisfy the equal time commuta- 
tion relations, 



< 

c 

3 




kA/n 

FIG. 3: Typical dispersion relation for coupled microres- 
onator system as depicted in Fig. 1 (solid line). For compar- 
ison, the dispersion relation of the system in the limit of no 
coupling (dotted line) is also shown. The resonance frequency 
of the cavity is given by LOonh/nc = 124.58 



g± OM) ,9± (z',t) = H z ~ z ') 

[&(*,*),&+(«',*)] = HZ-Z'): 



(18) 



with all other commutation relations vanishing. The 
function 5(z — z') is an effective delta function such that 
J-oo f( z )H z ~ z ')dz = f(z') when the function f(z) has 
its wavenumber restricted to the first Brillouin zone of 
the system. In terms of the effective fields, the Hamilto- 
nian in l|16|l becomes 



H = Tiuj b J dzg + g\ 



+ i 



- (-l) q hZ J dz(tfg-+c.c.) 





dg\ 




az 


9+ dz 




dg j 

-dz~ 9 - 


t dg- 
dz 




J dz (tfg. 


f. +C.C.) 





+ hujb I dzg^g_ — i — / dz — — g_ — gl-r — 



(19) 
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FIG. 4: Dispersion relation obtained using the transfer ma- 
trix technique (solid line) and the Hamiltonian in (flfifl (cir- 
cles), (a) The real part of the dispersion relation, (b) The 
imaginary part of the wavenumber for frequencies within the 
Bragg gap (c) The imaginary part of the wavenumber for fre- 
quency within the resonator gap. 



where uJb denotes the Bragg frequency centered at the 
Brillouin zone center and closest to wrjii. Using the 
Heisenberg equations of motion for the effective fields, 
we obtain the coupled equations 



d c d 
dt n dz 



' d_ 

dt 



cd_ 

n dz 



g+(z,t) = -iu b g + (z,t) + iEb(z,t) , 
g-(z,t) = -iuj b g-(z,t) + i(-l) q Zb(z,t), 



—b(z,t) = -iw b(z,t) + iZg + (z,t) + i (-1) 9 ^#_ (z,t) . 



(20) 



One can obtain the dispersion relation directly from 12U|) 
by assuming that each field is a plane wave e tkz - luJt j with 
k restricted to the first Brillouin zone. The results are 
equivalent to those in Fig. 4, obtained by diagonalizing 

Although the CME (|20|) were derived assuming an in- 
finite medium, they can be used to describe a struc- 
ture where the coupling constant S varies slowly over 
a distance on the order of the spacing between the res- 
onators. A multiple scale analysis can be used to iden- 
tify this limit and corrections to it. A more striking in- 
homogeneous structure is one beginning with a region 



(a) 100 




124.0 124.1 



t r 

124.2 124,3 



(t>) 100 




FIG. 5: Transmission spectrum for finite structure that con- 
tains 30 cavities, using parameters depicted in Table I. (a) 
Solid line represents the transmission spectrum obtained us- 
ing coupled mode equations and circles represents the trans- 
mission spectrum obtained using transfer matrix, (b) Trans- 
mission spectrum in the vicinity of the resonator gap using 
coupled mode equations (solid line) and transfer matrix (solid 
line with circles). 
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where there are no resonators, followed by a length L 
over which resonators are placed with an equal spacing 
and equal coupling to the channel(s), followed by a region 
where again there are no resonators. A simple model for 
such a region would be to use the equations 120f) . but 
replacing 3 with a position dependent coupling constant 



[8 (z) — 9 (z — L)] H, where 9 is the usual step function. 
It can be easily seen that this model formally violates 
our assumptions. Consider, for example, fields with a 
stationary time dependence, so g+(z,t) — g + (z) e~ luJt , 
and similarily for all other fields. Then the first equation 
gives 



- iug+ (z) + 



, dz 



1+ (z) = -iu;bg+ (z) +i[0 (z) — 8 (z — L)] Eb (z, t) , 



(21) 



r 



where in fact the factor [8 (z) — 8 (z — L)] could be omit- 
ted, since the third of (|20[1 together with the position 
dependent coupling constant guarantees that b {z) will 
only be nonzero in the region between z — and z = L. 
Note however that at z = and z — L the equation 121|) 
leads to a discontinuous dg + (z) /dz if it is assumed that 
g+(z) is everywhere continuous. This violates, of course, 
the assumption that fields such as g+ (z, t) are of the form 

03- 

Despite such a formal violation of our assumptions, 
this simple model in fact gives a good description of the 
optical response of a finite structure. To see this, consider 
first the fields g±(z,t) within the structure. It is clear 
from (|20|l that for a supposed frequency 57 there are two 
Bloch wavenumbers, which equivalently follow from QB) ; 
they are given by k{ui) — ±k, where 



n (A A!-S2) 2 -S4 



A 2 



(22) 



where g+, g_ are independent and q = um/c. For z < 
we denote the constants by g+ and <?^, and for z > L 
we denote them by g^ and g>. Now we consider the 
boundary condition at z — L, and note that since no field 
is incident from z > L, we have g^ = 0; an incident field 
is specified by g>f . Our independent unknowns are then 

g5,.g^, and the constants g+ and gq 2 that specify the 
field in the structure. We solve for these four unknowns 
by requiring the continuity of g±(z) at z — and z — 
L. The resulting transmittance of the structure can be 
written as 



with 



2» = 



> p iqL 



9+e 



9? 



,(2) 



(24) 



In the equation above Ao = (Q — uiq) is the detuning 
from the resonance frequency and Ai = (Q — uj(,) is the 
detuning from the Bragg frequency that lies closest to 
ujq. As a result, one can write the forward and backward 
propagating effective fields, g± (z,t), as 



(l) e 
9+ = ~ 



9+ = 



(2) _ e 



kA V S 2 



2 

ikL 



9le iqL . 



9± (z,t) 
9± ( z ) 



g± (z) e 



(2) - 

9± e 



■ikz 



(23) 



Once gv' and are set, g^ and g ( _^' are determined 
by the dispersion relation, or equivalently l|2(J|) . Hence 
there are only two independent constants. Outside the 
structure (3 = 0) there are also two independent con- 
stants in each of the regions z < and z > L, but the 
solution of H20|) is simpler. There it takes the form 



9±{z,t) 
9+(z) 
<?-(*) 



9±{z)e 

iqz 



= 9+e 
= 9-e~ 



In Fig. 5 we compare the transmission spectrum 
of a two channel micro-ring resonator structure with 
30 cavities, calculated both using the transfer matrix 
technique, 7 and using the coupled mode equation result 
eqn. 24. Again we adopt the parameters of Table I. Gen- 
erally there is good qualitative agreement, with the main 
features of the spectrum well described by the coupled 
mode equation result (24), although as noted above it is 
being applied beyond its strict range of applicability. An 
extension of this approach leads to the use of the CME 
(20) to treat a finite structure where the coupling con- 
stant S varies from one resonator to the next. To describe 
this we simply allow S in (20) to adopt a z-dependence, 
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c d 



di 



d_ 

di 



g+(z,t) = -ico b g+ (z, t) + iS (z) b (z, t) , 

-iuj b g^ (z, t) + i S (z) b (z, t) , 

b(z,t) = -ioj b(z,t) + iZ(z)g + {z,t)+i{-l) q Z{z)g-{z,t). 



d c d . 
at noz 



(25) 



r 



(a) lao- 




(b) ioo 




(oriA/nc 

FIG. 6: Transmission spectrum for short, finite, apodized 
structure with 5 unit cells, (a) Solid line represents the trans- 
mission spectrum obtained using transfer matrix and squares 
represent transmission spectrum obtained using coupled mode 
equations, (b) Transmission spectrum in the vicinity of the 
resonator gap using transfer matrix (solid line) and coupled 
mode equations (solid line with squares). 

In Fig. 6 we plot the transmission spectrum for 
a 5 cavity structure apodized such that the cavities 
(from left to right) are characterized by coupling con- 
stants (<7i,...cr 5 ) = (0.993,0.986,0.98,0.986,0.993), 
corresponding to (SiAn/7rc, ..^An/itc) = 
(0.0208,0.0287,0.0351,0.0287,0.0208). The transfer 



matrix results is presented, as well as a very simple 
application of the CME (25) using a piecewise uniform 
function to represent 3, where in the n th unit well we 
set S = S n . Again there is good qualitative agreement, 
although the CME are being applied beyond their strict 
range of applicability. Besides the difference between 
the CME and transfer matrix results with respect to 
the Fabry-Perot type oscillations, as seen in Fig. 5, 
here the CME solution also consistently overestimates 
the transmission on the high-frequency side of the stop 
gap. This can be traced back to the effects on the band 
curvature induced by the next highest Bragg gap, which 
are implicitly included in the transfer matrix solution 
but not in the CME calculation. 

Finally, we note that while at least three coupled mode 
equations are necessary to describe the kind of struc- 
tures we consider here if we deal with both their space 
and time dependence, if we instead restrict ourselves to a 
stationary time dependence, g± (z, t) = g± (z) e~ lult and 
b(z,t) = b(z)e~ ZUJt , then in fact we can eliminate the 
variable b (z, t) and construct coupled mode equations 
involving only g + {z,t) and <?_ {z,t). They are 



9 I \ 
d-z 9+ {Z) 

9 I \ 
oz 



w (uj) g+ (z) + i (-1) 9 n (u) g- (z) , 
-iv (u) g- (z) - i (-1) 9 ii (u>) g + (z) , 



(26) 



where 



^2 



(U}Q - Uj) 

n2 



(uj b - lu) 



-7—— — r, (27) 

C (CJQ - UJ) 

These equations are valid for u) ^ u)$. It is well-known 
that a photonic band gap opens in the dispersion relation 
described by these equations when > |^(oj)|, 21 

and that the width of the gap is larger for larger values 
of |/z(a->)|. Consequently we see from these equations an 
analytic confirmation of features that our dispersion re- 
lation display. Within our three mode model, one edge of 
the resonator gap occurs at ui — > ujq (in which case v and 
(i both diverge equally quickly and are hence equal in the 
limit as to approaches u>o), and one edge of the Bragg gap 
occurs at to — * u> b , because then the second term in the 
expression for v (to) vanishes, and v (u> b ) = [i (u) b ). 
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V. CONCLUSION 

We have presented a phenomena-logical Hamiltonian 
description of light propagation in side-coupled res- 
onators. This formulation is appealing in its simplicity, 
since it captures the basic physics of the structures via 
a set of readily understandably parameters. The most 
interesting special case is perhaps where a resonator gap 
is close to a Bragg gap, and at frequencies close to these 
gaps a three mode model gives a good description of the 
dynamics of a periodic structure of resonators. Coupled 
mode equations based on these captures the dispersion 
relation even deep within the gaps, and a naive extension 
of these equations to describe finite structures, although 
not within the strict range of applicability of the model, 
gives a good qualitative description. 

A hallmark of the kind of approach we have taken 
here is the connection of theoretically calculated or ex- 
perimentally observed parameters, such as the coupling 
coefficient a, to the parameters that appear in our phe- 
nomenological Hamiltonian. Such a strategy is particu- 
larly amenable to the description of quantum and non- 
linear optical effects. The Hamiltonian description leads 
to straightforward quantization, of course, and appropri- 
ate nonlinear terms can easily be added to the Hamilto- 
nian. In a previous study by Grimshaw et al. 22 , it was 
shown that three nonlinear coupled mode equations sup- 
port stationary solitary wave solution in the presence of 
Kerr nonlinearity. Numerical studies have indicated that 
soliton-like waves exist in resonator structures. In future 
work we plan to apply the approach we have detailed 
here to study such field excitations, where a Hamiltonian 
framework provides the ability to characterize conserved 
quantities in terms of the symmetries of the nonlinear 
field theory. 
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VII. APPENDIX 1 

In this appendix we use the Hamiltonian © to deter- 
mine the transmission properties of a single-cavity struc- 
ture. These transmission properties have been inten- 
sively studied using various methods such as finite dif- 
ference, time domain simulations—, and scattering ma- 
trix techniques 26 , and it is well-known that a Lorcntzian 
function gives an excellent approximation to the response 
of the structure. Here we show that our Hamiltonian also 
leads to a Lorentzian spectrum. To discuss transmission 
and reflection, we assume that there is a time-dependent 



source, u (t), coupled to the forward propagating modes 
at z s < Zq. We therefore modify the Hamiltonian J2J to 
include a source term: 



tng H source? 



(28) 



with 

H. source ^ 



a\u (t) e~ ikZs + a k u* (t) e 



ikz 



dk, 
(29) 

where e tkZs accounts for the fact that the light is gen- 
erated at z = z s . Using the Hamiltonian (|28|l and the 
commutation relations in the Heisenberg equations 
of motion (JSJ we find 

a k (t) = i& f b^e-^-^e-^dt' 
^ e -ia> h (t-t') e -ikz 3dt ^ 

/t 
-co 

-iu b(t) + i / tlkak{t)e l 



i u 

' — oo 



Cfe(i) 

db(t) 
dt 



- lkz °dk 



+ t(-i) q / ti- k c k {ty kz °dk. 



(30) 



where we have formally integrated the Heisenberg equa- 
tions for da k I dt and dc k / dt, so that both a k {t) and c k (t) 
are expressed entirely in terms of b (t) and u (t) . Using 
the expressions for a k (t) and c k (t) in the equation for 
db/dt, and expanding b (t) and u (t) in terms of Fourier 
components, 



b(t) 



u(t) = 



1 

2^ 
1 

27 



b {lu) e~ luJt duj, 

> 

u (w) e'^duj, 



we obtain 



-2tt<~/c 



u (to) e 1 



2im&,/c - iA 
where A = (lu — lj q + a (to)) and lu = luu/c, with 



a{uu) = 2 I p 

>0 



F 2 



LU 



dk 



(31) 



(32) 



(33) 



describing the small shift in the resonance frequency of 
the cavity due to the presence of the waveguide. To es- 
timate the effect of cx(lu), we assume takes a gaussian 
form in k space with a peak centered at k — luo. We take 
the width of the gaussian profile to be about l/im -1 , 
associated with a typical length over which the coupling 
between the waveguide and resonator is significant. Us- 
ing this approximate form for in the expression for 



11 



a{oj) and numerically evaluating the integral, we have 
verified that a(uS) is much smaller than the resonance 
frequency ujq for structures of interest. Note that in P2"|) 
we have switched our notation for wavenumber from k to 
u) = nuj/c — \k\ to stress that we are now considering the 
frequency response of the structure. To determine the 
transmission and reflection spectrum of the structure we 
define a set of effective fields 



/+(*,*) 



1 



dka k (t) e lkz 



/-CM) 



dkc k (t) e lkz 



(34) 



We then substitute the values 1)30(1 for a,k if) and Ck (t) in 
the effective fields l|34f) . and use the Fourier transforms 
(|31|l of b (t) and u (t) to simplify the integrals. We are 
specifically interested in the following two quantities 



lim /+ (z, t) 



-iA 



2wn4l/c-iA 



(Lu)e ik ( z - z ^e- luJt duj, 



z—> — oo 



1-Knillc 



(35) 



Note that in the absence of coupling we would have 



lim /+ (z, t) = - 

Z^OC C JQ 

lim /_ (z, t) = 



u (u>) e 



ik{z — z s ) —iujt 



du>, 



(36) 



The first (second) of the expressions in (|35H is the trans- 
mitted (reflected) field built as a superposition of the 
Fourier components of the source term, u(u>). We can 
therefore define the transmission and reflection coeffi- 
cients as 



t(w) 
r(w) 



-iA 



2irn£%/c- iA 
(-1) 



q 27T<g/c 



2tt<?/c- iA 



From these coefficients, it is clear that the cavity affects 
the transmission/reflection of the structure when the de- 
tuning, A, is on the order of 27rn£?/c. In the limit of very 
weak coupling - that is, when the value of 27m£?/c is ap- 
proximately constant over a frequency range centered at 
ujq and spanning several multiples of 27rn£?/c, then the 
transmission and reflection are well approximated by a 
Lorentzian lineshape 



t(u) 



-iA 



7 — iA 
r(w) * (-l) q 



7 — iA 



(37) 
(38) 



where 7 = 27rn£& /e. This condition yields 7 <^ cAk/2n\ 



for our assumed Ak ~ 27r/(l/xm) this gives the require- 
ment 7 -C 300ps _1 , which is met by typical values of 7 
(see equation ||5J and Table I). 



VIII. APPENDIX 2 

In this appendix we build the continuous collective op- 
erator bk (|12|) that applies for an infinite system of dis- 
crete resonators by first considering only excitations that 
are periodic over a length L — NA, and then passing to 
the limit N — > 00. In the periodic case there are still an 
infinite number of resonators, but only TV of the bi are 
independent. Assuming N is even, we can take them to 
be 



N 
~2 



1. 



N 

y 



N 
~2 



1. 



N 



(39) 



We denote this range by R. For an I outside Ri , we have 
bi = bi-pN where p is an integer such that 1—pN is within 
the range (|3^t . If we now introduce discrete wavevectors 
k m = 2-Km/L, where 



in 



N 
~2 



1. 



N 
~2 



N 
~2 



1 



TV 
2"' 



(40) 
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(that is, to £ i?) we can introduce Fourier amplitudes b r , 
according to 



b m = 



-ik m zi 



leR 



where zi = IK. We then find immediately that 



h = — L= V b 



mER 



and that 



while 



for example, so 



£&^ = £te 



£[^, 



or 



2tt 



£ 



= i 



(41) 



(42) 



(43) 



a form that we will presently find useful. 

We now consider letting N — > oo, with L — > oo such 
that A is fixed. Then the range i? approaches all the in- 
tegers from — ooto +oo, while k m become more closely 



spaced and approach a dense distribution of points rang- 
ing from — 7r/A to 7r/A; this is the first Brillouin zone, 
and we denote it by B.Z. In the usual way, then, we take 



2tt 



£ 



dk\ 



(44) 



B.Z. 



and. if we introduce bu such that 



bk, 



(45) 



where the k in b^ is first identified with k m but then 
allowed to vary continuously as N — ► oo, from (|43|l we 
have 



dk' 



'B.Z. 

and so we can identify 



h,bl] = 5{k-k'), 



for k and fe' within B.Z. In this limit, using I|44I45|) . we 
find 



B.Z. 



from (|42|l . where the integer Z now ranges from — oo to 
oo, and we recover l|12f> from l|41|) . 
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